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Employing mutual-inductance measurements we study the magnetic field dependence of com- 
plex AC susceptibility of artificially prepared highly ordered (periodic) two-dimensional Josephson 
junction arrays of unshunted Nb — AlOx — Nb junctions. The observed behavior can be explained 
assuming single-plaquette approximation of the overdamped model with an inhomogeneous critical 
current distribution within a single junction. 
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I. INTRODUCTION 

Despite the fact that Josephson Junction Arrays (JJA) have been actively studied for decades, they continue to 
contribute to the variety of intriguing and peculiar phenomena. To name just a few recent examples, it suffice to 
mention the so-called paramagnetic Meissner effect and related reentrant temperature behavior of AC susceptibility, 
observed both in artificially prepared JJA and granular superconductors (for recent reviews on the subject matter, see, 
e.g. [1-4] and further references therein). So far, most of the investigations have been done assuming an ideal (uniform) 
type of array. However, it is quite clear that, depending on the particular technology used for preparation of the array, 
any real array will inevitably possess some kind of non-uniformity, either global (related to a random distribution of 
junctions within array) or local (related to inhomogeneous distribution of critical current densities within junctions). 
For instance, recently a comparative study of the magnetic remanence exhibited by disordered (globally nonuniform) 
3D- JJA in response to an excitation with an AC magnetic field was presented [5] . The observed temperature behavior 
of the remanence curves for arrays fabricated from three different materials {Nb, Y Ba2Cu30'j and Lai.85<S'ro.i5Cu04) 
was found to follow the same universal law regardless of the origin of the superconducting electrodes of the junctions 
which form the array. 

In the present paper, through an experimental study of complex AC magnetic susceptibility x(2^; ^^c) of the 
periodic (globally uniform) 2D-JJA of unshunted Nb — AlOx — Nb junctions, we present evidence for existence of 
the local type non-uniformity in our arrays. Specifically, we found that in the mixed state region x{T, fiAc) can 
be rather well fitted by a single-plaquette approximation of the overdamped 2D- JJA model assuming a nonuniform 
(Lorentz-like) distribution of the critical current density within a single junction. 

II. EXPERIMENT 

Our samples consisted of 100 x 150 unshunted tunnel junctions. The unit cell had square geometry with lattice 
spacing a — A6^m and a junction area of 5 x The critical current density for the junctions forming the arrays 

was about 6Q0A/cm^ at 4:.2K, giving thus Ic = 150/iA for each junction. 

We used the screening method [6] in the reflection configuration to measure the complex AC susceptibility x — x'+x" 
of our 2D- JJA (for more details on the experimental technique and set-ups see [7-9] ). Fig.l shows the obtained 
experimental data for the complex AC susceptibility x(r, Hac) as a function of the excitation field Hac for a fixed 
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temperature below Tc- As is seen, below 50mOe (which corresponds to a Meissner-like regime with no regular flux 
present in the array) the susceptibility, as expected, practically does not depend on the applied magnetic field, while 
in the mixed state (above 50mOe) both x'(T, /i^c) and x"(T, /i^c) follow a quasi-exponential field behavior of the 
single junction Josephson supercurrent (see below). 



III. DISCUSSION 



To understand the observed behavior of the AC susceptibility, in principle one would need to analyze the flux 

dynamics in our ovcrdamped, unshunted JJA. However, given a well-defined (globally uniform) periodic structure 
of the array, to achieve our goal it is suflicient to study just a single unit cell (plaquette) of the array. (It is worth 
noting that the single-plaquette approximation proved successful in treating the temperature reentrance phenomena 
of AC susceptibility in ordered 2D- JJA [4,7,8] as well as magnetic remancncc in disordered 3D- JJA [5].) The unit 
cell is a loop containing four identical Josephson junctions. Since the inductance of each loop is £ = hqu ~ QApH 
and the critical current of each junction is Ic = 150/uA, for the mixed-state region (above bOmOe) we can safely 
neglect the self- field effects because in this region the inductance related flux ^L{t) = ^I{t) (here I{t) is the total 
current circulating in a single loop [10]) is always smaller than the external field induced fiux $ea:t(i) = Bac{t)S (here 
5 ~ is the projected area of a single loop, and Bac{t) = ixoHac cosut is an applied AC magnetic field). Besides, 
since the length L and the width w of each junction in our array is smaller than the Josephson penetration depth 
\j = \/ ^o/'^T^I^odjcO (where jco is the critical current density of the junction, <l>o is the magnetic flux quantum, and 
d = 2Ai -|- ^ is the size of the contact area with Xl{T) being the London penetration depth of the junction and ^ 
an insulator thickness), namely L w 5/iTO and Aj ~ 20/iTO (using jco = GOOA/cm^ and = 39nm for Nb 
at T = 4.2K), we can adopt the small-j unction approximation [10] for the gauge-invariant superconducting phase 
difference across the ith junction (by symmetry we assume that [7,8] 4>i = (1)2 = 4>3 = = ^i) 

27TBac{t)d. 

(t)i{x, t)=(l)Q-\ — ^ X (1) 

^0 

where (f)o is the initial phase difference. 

The net magnetization of the plaquette is M{t) = SIs{t) where the maximum supercurrent (corresponding to 
'Po = T^/^) through an inhomogeneous Josephson contact reads 

L w 

Is{t) = j dx j dyjc{x,y)cos(j)i{x,t) (2) 



For the explicit temperature dependence of the Josephson critical current density 



jco(T)=jco(0) 



A(T) 
A(0) 



tanh 



A(r) 



2kBT 



(3) 



we used the well-known [11] analytical approximation for the BCS gap parameter (valid for all temperatures), A(T) = 

A(0)tanh ^2.2^^^^ with A(0) = 1.76fcsTc. 

In general, the values of x'{T,hAc) and ^' {T,hAc) of the complex harmonic susceptibility are defined via the 
time-dependent magnetization of the plaquette as follows: 

X'(T, hAc) = -T — / d{ujt) cos{Lut)M{t) (4) 

T^flAC Jo 

and 

x"{T,hAc) = ^— [ d{ujt)sm{ujt)M{t) (5) 

T^rlAC Jo 

Using Eqs. (l)-(5) to simulate the magnetic field behavior of the observed AC susceptibility of the array, we found 
that the best fit through all the data points and for all temperatures is produced assuming the following nonuniform 
distribution of the critical current density within a single junction [10] 
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ic(.,.)=ico(T)(^)(^) (6) 

It is worthwhile to mention that in view of Eq.(2), in the mixed-state region the above distribution leads to approxi- 
mately exponential field dependence of the maximum supercurrent /s(T, hAc) — IsiT, 0)e~'^^^^'^° which is often used 
to describe critical-state behavior in type-II superconductors [12]. Given the temperature dependencies of the Lon- 
don penetration depth Xl{T) and the Josephson critical current density jco{T), we find ho{T) = ^o/2tthoXj{T)L ~ 
ft-o(0)(l — T/Tc)^/^ for the temperature dependence of the characteristic field near Tc- This explains the improvement 
of our fits (shown by solid lines in Fig.l) for high temperatures because with increasing the temperature the total 
flux distribution within a single junction becomes more regular which in turn validates the use of the small-junction 
approximation. 

In summary, we found clear experimental evidence for the influence of the junction nonuniformity on magnetic 
field penetration into the periodic 2D array of unshunted Josephson junctions. By using the well-known AC magnetic 
susceptibility technique, we have shown that in the mixed-state regime the AC field behavior of the artificially prepared 
array is reasonably well fitted by the single-plaquettc approximation of the overdamped model of 2D-JJA assuming 
inhomogeneous (Lorentz-like) critical current distribution within a single junction. 



ACKNOWLEDGMENTS 



We thank P. Barbara, C.J. Lobb, and R.S. Ncwrock for useful discussions. S.S. and F.M.A.M. gratefully acknowledge 
financial support from Brazilian Agency FAPESP under grant 2003/00296-5. 



[1] R.S. Newrock, C.J. Lobb, U. Geigenmuller and M. Octavio. Solid State Physics 54 (2000), p. 263. 
[2] P. Martinoli and C. Leeman. J. Low Temp. Phys. 118 (2000), p. 699. 

[3] S. Sergeenkov. In: A. Narlikar, Editor, Studies of High Temperature Superconductors, vol. 39, Nova Science, New York 
(2001), p.ll7. 

[4] F.M. Araujo-Moreira, P. Barbara, A.B. Cawthorne and C.J. Lobb. In: A. Narlikar, Editor, Studies of High Temperature 
Superconductors, vol. 43, Nova Science, New York (2002), p. 227. 

[5] S. Sergeenkov, W.A.C. Passos, P.N. Lisboa-Fillio and W.A. Ortiz. Phys. Lett. A 291 (2001), p. 311. 

[6] J.L. Jeanneret, G.A. Gavilano, A. Racine, Ch. Leemann and P. Martinoli. Appl. Phys. Lett. 55 (1989), p. 2336. 

[7] F.M. Araujo-Moreira, P. Barbara, A.B. Cawthorne and C.J. Lobb. Phys. Rev. Lett. 78 (1997), p. 4625. 

[8] P. Barbara, F.M. Araujo-Moreira, A.B. Cawthorne and C.J. Lobb. Phys. Rev. B 60 (1999), p. 7489. 

[9] W.A.C. Passos, P.N. Lisboa-Filho, R. Caparroz, C.C. de Faria, P.C Venturini, F.M. Araujo-Moreira, S. Sergeenkov and 
W.A. Ortiz. Physica C 354 (2000), p. 189. 
[10] A. Barone and G. Patcrno. Physics and Applications of the .Josephson Effect, Wiley, New York (1982). 
[11] R. Meservey and B.B. Schwartz. In: R.D. Parks, Editor, Superconductivity, vol.1, M. Dekker, New York (1969), p. 117. 
[12] F.M. Araujo-Moreira, O.F. de Lima and W.A. Ortiz. J. Appl. Phys. 80 (1996) p. 6. 



3 




10 100 



0,06- 


(b) 






0,05- 








0,04- 




/ ' 




0,03- 








0,02- 








0,01 - 









10 100 

h,, (mOe) 



0,06- 


(d) 






0,04- 








M 








0,02- 








0,00- 










10 100 



4 




FIG. 1. The dependence of the magnetic susceptibihties, x' (T, hAc) and x"{T,hAc), on AC magnetic field amplitude h^c 
for different temperatures: T = 4:.2K (a,b), T — 6K (c,d), and T = 87^ (e,f). Solid lines correspond to the fitting of the 
2D-JJA model with nonuniform critical current profile for a single junction (see the text). 
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